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I.  INTRODUCTION 

In  the  work  described  in  this  report,  the  phenomenon  of  target  radar 
cross  section  dependence  on  transmitted  frequency  was  investigated  as  a  can¬ 
didate  technique  for  target  classification*  Very  simply,  the  phenomenon 
results  from  the  fact  that  a  target  consists  of  several  individual  reflectors 
which  constructively  or  destructively  interfere,  depending  on  their  range 
separation  and  on  the  tranmitted  frequency*  Thus,  over  a  given  frequency 
diverse  bandwidth,  the  received  signal  amplitude  will  fluctuate  or  scin¬ 
tillate,  and  these  scintillations  can  be  analyzed  to  determine  target 
characteristics. 

Frequency  diversity  has  long  been  investigated  as  a  method  of  improving 
target  detection1 -4,  tracking  accuracy  5-7 f  an(j  range  resolution.®  More 
recently,  frequency  diverse  target  scintillation  has  been  studied  as  a  tech¬ 
nique  for  acquisition  and  classification  of  stationary  targets. 9-13  Also,  fre¬ 
quency  diversity,  combined  with  polarization  diversity,  has  been  investigated 
both  experimentally  and  theoretically;  in  these  techniques,  amplitudes  and 
relative  phases  are  measured  to  obtain  stationary  target  signatures.1^-1® 

The  purpose  of  the  present  experimental  effort  was  to  design  and  build  a 
94  GHz  radar  and  to  evaluate  the  utility  of  the  radar  for  detecting  frequency 
diverse  radar  cross  section  effects.  A  few  additional  measurements  were  made 
with  a  140  GHz  radar,  and  these  data  are  also  presented  and  discussed. 

In  Section  II,  a  brief  theoretical  analysis  of  expected  results  is 
outlined.  Descriptions  of  the  two  radar  systems  used  in  this  experiment  are 
given  in  Section  III,  and  the  signal  processing  and  measurement  techniques  are 
also  discussed. 

Experimental  data  taken  for  corner  reflectors  at  various  range  spacings 
are  presented  in  Section  IV.  The  frequency  calibration  method  used  1b 
described,  and  range  separation  values  obtained  with  radar  results  are  com¬ 
pared  with  measured  values. 

For  the  94  GHz  data,  a  computer  simulation  program  was  developed  for  com¬ 
parison.  A  description  of  the  program  is  given  in  Section  V,  while,  in 
Section  VI,  the  agreement  and  discrepancies  between  the  theoretical  and 
experimental  results  are  discussed. 

Finally,  in  Section  VII,  conclusions  from  the  simulated  and  experimental 
results  are  drawn  with  respect  to  radar  characteristics,  and  recommendations 
for  further  research  are  given. 

II.  MATHEMATICAL  ANALYSIS 

Consider  a  target  consisting  of  a  reflectors  with  radar  cross  sections 
(RCS)  c; ,  02,  03,  .  .  .  on.  Assume  that  the  target  is  stationary  and  that  all 
reflectors  are  equally  illuminated  by  the  radar.  The  target  is  in  the  far- 
fleld  of  the  radar,  and  vice  versa.  Assume  also  that  the  typical  linear 
dimension  of  the  reflectors  is  much  larger  than  the  transmitted  wavelength 
(optical  region).  In  general,  the  RCS  of  the  individual  reflectors  will  vary 
as  the  frequency  squared,  as  for  flat  plates  and  corner  reflectors.  However, 
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the  frequency  diverse  bandwidth  is  assumed  to  be  sufficiently  small  that  these 
changes  in  RCS  are  negligible. 

Let  be  the  range  to  the  ith  reflector,  and  define 

®i  -  -  Rij  .  (D 

that  is,  R34  is  the  range  separation  between  the  ith  and  the  jth  reflector. 

It  can  be  shown  that,  for  a  continuous  ware  radar,  the  square-law  detected 
signal  voltage,  V(f),  is  given  by^ 


V(f)  -  A  {c;  +  02  +  °3  + . +  °n 


Here  A  is  a  parameter  depending  on  target  range,  transmitted  power,  and 
characteristics  of  the  radar;  c  is  the  speed  of  light,  and  f  is  the 
transmitted  frequency. 

If  the  voltage  signal  is  measured  as  a  function  of  frequency,  there  is  a 
constant  "dc"  level  and  a  superposed  "ac"  waveform  which  is  periodic  with 
frequency. 

The  result  is  analogous  to  a  time  domain  periodic  wave  form;  here  time  is 
analogous  to  frequency,  and  frequency  is  analogous  to  the  quantitiy 


Rote  that  this  "frequency"  f^j  is  the  time  it  takes  electromagnetic  energy  to 
make  a  round  trip  between  the  ith  and  Jth  reflectors. 


Equation  (3)  shows  that,  mathematically,  the  Fourier  transform  of  V(f)  is 
already  performed.  The  transform  of  V(f)  should  show  a  spike  of  amplitude 

Aij  “  2A^Oicrj  (5) 

at  a  "frequency"  of  f^j  (Equation  (4)). 

Each  pair  of  reflectors  contributes  a  spike  to  the  Fourier  transform, 
although  two  pairs  with  equal  range  separation  contribute  to  the  same  spike 
and,  thus,  cannot  be  distinguished.  However,  a  high  degree  of  structural 
periodicy  in  the  target,  i.e.,  many  reflectors  with  equal  spacing,  would 
result  in  a  strong  spike  at  the  characteristic  range  separation.  The  number 
of  contributions  to  the  Fourier  transform  spectrum  increases  in  an  arithmetic 
series  with  the  number  of  reflectors;  for  n  reflectors,  there  are  1^?n(n-1 ) 
contributing  pairs. 

All  the  experimental  data  taken  were  for  two  or  three  reflectors.  For 
the  two  reflector  case,  the  V(f)  (or  "interference  pattern")  is  given  by 


V(f)  -  A{o;  ♦  o2  ♦  cos[2u(^)  f]}  (6) 

i.e.,  a  single  sinusoidal  oscillation.  For  three  reflectors, 

V(f)  -  A{o;  +  02  +  °3 

+  2^0; op  cos[2n^l2j  f] 

+  2^0;  Cj  COS  [2ir^2R^  fj 

+  ajojjoj  cos[2n^22^  f]}  (7) 

so  that  V(f)  in  this  case  is  a  superposition  of  three  sinusoidal  oscillations. 

In  terms  of  these  oscillations,  one  complete  cycle  will  be  observed  when 
the  cosine  argument  changes  by  2n,  or  when  the  change  in  frequency,  Af,  is 
such  that 


if  ■  55^.  (8) 

The  frequency  diverse  bandwidth  (Af)  required  to  see  one  complete  oscillation 
in  V(f)  is  calculated  in  Table  1  for  various  reflector  separations.  The 
values  were  chosen  to  aid  in  understanding  the  experimental  data. 
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TABLE  1.  FREQUENCY  DIVERSE  BANDWIDTHS 


Reflector  Range 
Separation  Rij 
(Meters) 

Frequency  Diverse  Bandwidth 
Required  for  One  Cycle  in  V(f) 
(MHz) 

10 

15.0 

6 

25.0 

5 

30.0 

4 

37.5 

3 

50 

2 

75 

1 

150 

0.5  (50  cm) 

300 

0.1  (10  cm) 

1.5  GHz  1 

0.01  (1  cm) 

15  GHz 

III.  DESCRIPTION  OF  RADAR  SYSTEMS 

A.  94  GHz  Radar  System 

The  94  Ghz  radar  system  is  depicted  in  Figure  1 .  The  millimeter 
source  is  a  continuous  wave  (cw)  klystron  with  a  nominal  power  of  300  mW  and  a 
nominal  electronic  tuning  3  dB  bandidth  of  200  MHz.  The  klystron  reflector 
voltage  is  sawtooth  wave  modulated  at  a  1  KHz  repetition  rate  to  sweep  the 
klystron  over  its  entire  tuning  range.  Some  of  the  output  is  coupled  to  a 
thermistor  for  monitoring  the  power.  Also,  some  of  the  power  is  coupled 
through  a  cavity  wavemeter  (for  ballpark  frequency  measurements)  to  a  square 
law  detector  (for  monitoring  the  transmitted  power  waveform).  The  power  at 
the  transmit  antenna  for  this  experiment  is  estimated  to  be  about  200  mW.  A 
linear  vertically  polarized  wave  was  transmitted. 

The  receiver  consists  of  an  antenna,  a  square  law  detector  with  a  sen¬ 
sitivity  of  about  600  mV/mW,  and  a  low  noise  preamplifier  used  at  a  gain  of 
ten  and  an  input  impedance  of  1  MegOhm.  The  receive  signals  were  monitored 
with  a  low  bandwidth  (500  KHz)  high  sensitivity  (lOO*V/division)  oscilloscope 
triggered  by  a  1  KHz  sync  pulse  from  the  waveform  generator.  As  shown  in 
Figure  1 ,  the  received  waveform  consisted  of  the  broad  "interference  pattern" 
(received  signal  vs  klystron  reflector  voltage)  during  the  sawtooth  rise,  and 
a  compressed  version  during  the  sawtooth  fall. 

This  signal  was  applied  to  a  boxcar  averager/gated  integrator,  also 
triggered  at  1  KHz  by  the  waveform  generator.  The  boxcar  averager  provides  a 
variable  size  aperature  which  may  be  positioned  anywhere  during  the  received 
interference  pattern.  The  gated  integrator  is  operated  in  the  exponential 
averaging  mode,  which  provides  a  weighted  average  output;  that  is,  the  most 
recent  pulses  have  more  influence  than  early  pulses.  The  output  of  the  boxcar 
averager  can  be  considered  to  be  the  average  of  2  TC/AD  repetitions  of  the 
signal  falling  within  the  aperture,  where  TC  is  the  selected  integration  time 
constant  and  AD  is  the  selected  aperture  duration.  For  this  experiment, 

TOI.O  s  and  AD-10  s,  so  there  was  essentially  no  pulse  integration. 
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Boxcar  Averaged  Output 


The  boxcar  averager  provides  for  the  automatic  scan  of  the  aperture  bet¬ 
ween  selected  limits  and  also  provides  a  scan  voltage  proportional  to  the 
aperture  position.  Thus,  data  were  taken  by  sweeping  the  aperature  once 
across  the  broad  interference  pattern  and  plotting  boxcar  integrator  output  vs 
scan  voltage  (or,  equivalently,  reflector  voltage)  on  an  X-Y  recorder.  The 
sweep  time  was  about  18  s,  and  a  boxcar  integrator  signal  time  constant  of  0.1 
s  (low  pass  filter)  provided  some  smoothing  of  the  interference  pattern. 

B.  140  GHz  Radar  System 

The  characteristics  of  the  140  GHz  radar  have  been  documented  exten¬ 
sively  and  specifications  are  presented  in  Tables  2-4.^  The  frequency  of  the 
radar  is  manually  adjustable  over  a  600  MHz  bandwidth.  Data  were  taken  by 
measuring  the  pulse  height  of  the  target  response  on  an  oscilloscope  (A-scope) 
as  a  function  of  frequency.  Only  the  observed  maxima  and  minima  points  were 
recorded. 


TABLE  2.  TRANSMITTER  PERFORMANCE  SPECIFICATIONS. 


f  TRANSMITTER  -  Coherent/Pulse  Compression  | 

Operating  Frequency 

142.02  GHz 

Transmitted  Peak  Power 

10  mW  (10  dBm)  j 

Equivalent  Peak  Power 

600  mW  (28  dBm) 

Transmitted  Pulse  Width 

3  usee 

Compressed  Pulse  Width 

50  nsec 

Chirp  Bandwidth 

20  MHz 

Pulse  Compression  Ratio 

60 

PRF 

30  kHz 

Pulse  Repetition  Interval 

33*3  usee 

Duty  Cycle 

3% 

Average  Transmitted  Power 

0.9  mW 

Frequency  Diversity  Bandwidth 

660  MHz 

TABLE  3.  RECEIVER  PERFORMANCE  SPECIFICATIONS. 


RECEIVER  -  Pulse-Doppler  with 

Pulse  Compression 

First  LO  Frequency 

140i67  qHz  - 

First  IF  Frequency 

1395  MHz  1 

Second  LO  Frequency 

1275  MHz 

Second  IF  Frequency 

120  MHz 

Noise  Figure  (SSB)  of  First  Mixer 

13  dB 

First  Mixer  Preamplifier  Gain 

26  dB  minimum 

Circulator  Frequency  Range 

142+1  GHz 

Circulator  Isolation 

18  dB  minimum 

Circulator  Insertion  Loss 

1.6  dB  maxi  mum 

Circulator  VSWR 

1 .40  maximum 

Unambiguous  Range 

~5  km 

Receiver  Bandwidth 

24  MHz 

SAW  Line  Bandwidth 

25  MHz 

Range  Cell  (3-usec  pulse) 

1500  ft 

Range  Cell  (50-nsec  pulse) 

25  ft 
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TABLE  4.  TRACKING/GUIDANCE  ANTENNA  SPECIFICATIONS. 


Type 

Cassegrain,  Conscan 

Operating  Frequency 

142  GHz 

Diameter  (1) 

60  cm 

3-dB  Beamwidth  (0g  -  A/A) 

0*23°  ■  4*0  mrad 
(0.23°  -  4*0  mrad) 

Gain  (GD  "  4*  A/x2) 

56.2  dB  (59  dB) 

(-*■  efficiency  -  505?) 

Sidelobe  Level 

-16.3  dB 

Crossover  Points 

3  dB  «■+  2.0  mrad  offset 

1.5  dB  «-+  1.1  mrad  offset 

VSVR 

1.15 

Polarization 

Linear,  Vertical  ! 

Conical  Scan  Rate 

30  to  90  Hz 

Angle  of  First  Sidelobe 

(-6  mrad) 

Far  Field  (R  >  2D2/X) 

(R  >  340  m) 

IV.  EXPERIMENTAL  RESULTS 

A.  94  GHz  Radar  Data 

The  measurements  made  with  the  94  GHz  radar  were  taken  in  the  High 
Bay  area  of  Building  5400,  Redstone  Arsenal,  Alabama.  The  target  range  was 
-80  meters,  satisfying  far- field  conditions  for  both  radar  and  target.  The 
transmit  and  receive  antennas  were  side-by-side,  and  their  centers  were 
approximately  1  meter  above  the  concrete  floor.  Antenna  beamwidths  are 
0.4°  and  0.6°,  so  that  the  antenna  main  lobes  did  not  intersect  the  floor  at 
the  target  (no  main  lobe  multipath).  The  3  dB  beamwidths  of  the  transmit  and 
receive  antennas  overlapped  each  other  by  approximately  30  cm  at  the  target. 

The  first  data  taken  were  a  series  for  the  two  corner  reflectors 
separated  by  6,  5,  4,  3,  2,  and  1  meters.  These  results  are  shown  in  Figures 
2-7.  The  corner  reflectors  were  triangular  trihedrals  with  RCS  of  200  sm  and 
91  sm.  The  larger  corner  reflector  was  behind  and  above  the  smaller  corner 
reflector.  The  scales  shown  on  the  data  are  relative  amplitudes 
(uncalibrated);  the  klystron  reflector  voltage  range  is  about  100  V,  while  +he 
video  output  (boxcar  integrator  output)  is  on  the  order  of  millivolts.  These 
data  show  the  expected  result— a  single  oscillation  modified  by  the 
transmitter  power  curve.  As  a  measure  of  the  transmitter  (and  receiver)  power 
curve,  data  were  taken  for  the  transmit  waveform  monitor  output  (Figure  8) 
and  for  each  corner  reflector  separately  (Figure  9  and  10).  The  power  curve 
of  Figure  8  is  similar  to  typical  reflex  klystron  power  output  characteristics 
(Figure  II).20  Using  the  boxcar  aperature  in  conjunction  with  the  cavity 
wavemeter,  the  total  operating  range  was  determined  to  be  approximately  from 
93*96  Ghz  to  94.11  Ghz  (+  0.01  GHz),  or  a  total  of  about  150  MHz.  The  center 
frequency  (maximum  power  point)  was  94*05  GHz.  The  3  dB  electronic  tuning 
bandwidth  is  around  100  MHz,  which  is  about  half  that  claimed  by  the  manufac¬ 
turer. 


The  klystron  frequency  is  clearly  not  a  linear  function  of  refleotor 
voltage,  so  the  data  obtained  for  2  reflectors  separated  by  6  meters  was  used 
to  calibrate  frequency  change  vs  reflector  voltage  change.  These  calibration 
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Figure  4.  Two  reflectors  separated  by  4  meters. 
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Klystron  Reflector  Voltage  Change 
6.  Two  reflectors  separated  by  2  meters 
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Figure  11.  Power  output  characteristics  of  a  reflex  klystron  (after  Liao) 
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data  were  then  used  to  analyze  the  results  for  other  separations  and,  in  par¬ 
ticular,  compare  the  range  separation  calculated  from  the  radar  data  to  the 
actual  range  separation  as  determined  with  a  tape  measure. 


Figure  12  indicates  how  the  calibration  was  made.  The  position  of  each 
maxima  and  minima,  indicated  by  a  vertical  line,  was  estimated  to  the  nearest 
one  half  linear  unit  by  bisecting  the  given  half  cycle  at  its  half  height. 
These  data  are  given  in  Table  5.  From  Table  1,  it  is  seen  that,  for  a  6  meter 
separation,  the  maxima  are  25  MHz  apart,  so  that  the  extrema  are  12.5  MHz 
apart.  The  first  maximum  was  arbitrarily  assigned  a  value  of  25  MHz.  The 
data  of  Table  5  are  plotted  in  Figure  13;  typical  frequency  characteristics 
shown  in  Figure  14  are  similar. 20 

TABLE  5*  FREQUENCY  CHANGE  CALIBRATION  DATA 


Extrema 

Klystron  Reflector 
Voltage  Change 
(linear  units) 

Change  In 
Frequency 
(MHz) 

Maximum  #1 

11.0 

25-0* 

Minimum  #1 

18.0 

37.5 

Maximum  #2 

27-0 

50.0 

Minimum  #2 

34.5 

62.5 

Maximum  #3 

42.5 

75-0 

Minimum  #3 

51-5 

87.5 

Maximum  #4 

60.0 

100.0 

Minimum  #4 

67.0 

112.5 

Maximum  #5 

72.5 

125.0 

Minimum  #5 

79-0 

137.5 

Maximum  #6 

83.0 

150.0 

Minimum  #6 

87-5 

162.5 

♦Arbitrarily  Assigned 

For  the  data  in  Figures  3-7,  the  number  of  cycles  observed  and  the  reflector 
voltage  change  for  these  cycles  were  determined.  The  corresponding  frequency 
change  was  then  found  from  Figure  13 >  and  the  distance  between  reflectors  was 
calculated  from  Equation  (8)  in  the  form 


^•2M.^Af  -  (#  of  cycles)  (9) 

or 

AR  -  c(#  of  cycles)  (10) 

2AF 

where  Af  is  now  the  frequency  change  over  the  number  of  cycles  observed.  The 
results  of  the  radar  data  are  compared  with  the  actual  separations  in  Table  6. 
The  actual  separations  were  determined  with  a  tape  measure  between  the  front 
edges  of  the  corner  reflectors*  Because  the  phase  centers  of  the  reflectors 


Figure  1A.  Frequency  characteristics  of  a  reflex 
klystron  (after  Liao) . 
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are  not  precisely  known,  the  estimated  error  for  this  measurement  is  *_  10  cm. 
The  agreement  is  quite  good,  better  than  the  estimated  error,  except  for  the 
one  meter  separation  case,  where  the  power  curve  severely  distorts  the  inter¬ 
ference  pattern. 

TABLE  6.  REFLECTOR  SEPARATION  DETERMINED  FROM  RADAR  DATA 


Measured 

Separation 

Number  of 
Cycles  Observed 

Frequency 

Change 

(MHz) 

Calculated 

Separation 

(meters) 

Percent 

Error 

5-0 

41/£ 

136 

4-96 

oTejT 

4-0 

31/2 

132 

3.98 

0.5% 

3-0 

21/2 

125.5 

2.99 

0.3% 

2.0 

H/2 

110.5 

2.04 

2.0% 

1 .0 

_ 1/2 _ 

51 

1.47 

47.0 

An  explanation  of  the  fine  structure  (small  oscillations)  in  some  of  the 
data,  particularly  Figures  5,  6,  7,  9,  and  10,  is  required.  The  large  (~30  ft 
by  35  ft) ,  flexible  metal  door  of  the  High  Bay  was  approximately  20  meters 
behind  the  corner  reflectors.  The  radar  boresight  was  at  an  oblique  angle  to 
the  door,  so  that  under  normal  conditions  no  energy  was  reflected  from  the 
door  back  to  the  radar.  However,  when  the  door  was  buffeted  by  the  wind, 
bringing  it  more  perpendicular  to  boresight,  its  effect  could  clearly  be  seen 
in  the  radar  signal.  The  period  of  the  oscillations  is  consistent  with  the 
distance  between  the  door  and  target. 

Two  interference  patterns  were  recorded  for  a  target  consisting  of  three 
corner  reflectors.  A  50  am  reflector  was  placed  between  the  two  initially 
used.  For  the  first  pattern  (Figure  15)*  the  spacing  was  4  meters  between  the 
first  and  second  and  2  meters  between  the  second  and  third;  for  the  second 
pattern  (Figure  16),  the  spacing  was  3  meters  and  2  meters.  These  patterns, 
clearly  more  complex  than  those  for  two  reflectors,  are  analyzed  and  compared 
with  simulated  data  in  Section  6. 

B.  140  GHz  Radar  Data 

The  measurements  made  with  the  140  GHz  radar  were  taken  at  the  Road 
Test  Course  Area,  Redstone  Arsenal,  Alabama.  The  target  range  was  about  480 
meters,  and  two  triangular  trihedral  corner  reflectors  of  approximately  5,000 
am  and  6,000  sm  were  used.  Only  the  extrema  were  measured  on  an  oscilloscope 
while  manually  changing  the  transmitted  frequency.  The  frequency  calibration 
method  is  described  elsewhere. 19 

The  interference  patterns  for  two  different  reflector  separations  are 
shown  in  Figure  17  and  Figure  18.  Calculations  of  the  radar  range  separation 
of  the  corner  reflectors  using  the  radar  data  gave  1.61  meters  and  2.85 
meters;  both  values  agree  with  the  measued  values  within  about  5%.  Clearly, 
the  data  suffer  from  nonlinearities  in  the  receiver  or  transmitter.  This 
effect  is  further  demonstrated  in  the  response  for  a  single  corner  reflector 
(Figure  19),  taken  at  a  range  of  400  meters.  It  is  expected  that  such  extreme 
nonlinearities  will  make  analysis  of  results  for  more  than  two  reflectors 
essentially  impossible. 
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.  Three  reflectors  separated  by  3m,  2m. 
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V.  COMPUTER  SIMULATION 


A  computer  simulation  was  developed  to  compare  theoretical  versus 
measured  results.  Figure  20  details  the  sequence  of  actions  in  the  program 
which  develops  the  theoretical  response  of  the  square  law  detector  as  the 
voltage  of  the  source  is  swept.  Since  the  frequency  output  of  the  source  is 
not  linear  with  the  applied  slewing  (reflector)  voltage,  the  plots  are  ampli¬ 
tude  versus  slewing  voltage.  Figure  21  provides  the  sequence  of  actions 
resulting  in  the  discrete  Fourier  transform  (DFT)  of  the  theoretical  result. 
The  FFT  algorithm  used  to  compute  the  DFT  has  been  modified  to  compute  the  DFT 
on  a  computer  (HP  #9830)  which  does  not  have  a  COMPLEX  statement.  The  modifi¬ 
cation  divides  all  operations  and  coefficients  into  real  and  imaginary  parts. 

The  HP  #9830  uses  a  type  of  BASIC  (peculiar  to  HP  machines).  The  actual 
programs  are  shown  in  Appendices  A  and  B.  A  simple  way  of  inputing  the  N 
square  law  responses  is  to  use  the  LINK  command  in  the  appropriate  spot  in  the 
square  law  detector  response  program  (after  statement  number  620).  This  is 
necessary  for  a  small  machine  such  as  the  HP  #9830,  since  the  programs  are  on 
cassette  tape  and  no  statement  equivalent  to  the  FORTRAN  CALL  statement  exists 
in  the  HP  BASIC  language  used.  The  program  in  Appendix  C  has  the  LINK  com¬ 
mand,  directing  the  program  to  the  FFT  routine  stored  as  program  2  on  the  l 
cassette  tape.  In  order  to  use  this  feature,  the  number  of  frequency  steps 
must  be  a  power  of  2.  The  maximum  number  of  steps  allowed  by  the  structuring 
of  the  machine  is  256. 

VI.  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 

For  the  case  of  two  reflectors,  Figures  22  through  29  show  a  comparison 
of  simulation  (dashed  curves)  and  experimental  results  (solid  curves).  The 
simulation  results  shown  include  weighting  by  the  power  spectrum  shown  in 
Figure  8  and  the  nonlinear  frequency  dependence  on  reflector  voltage.  The 
frequency  of  modulation  produced  by  the  simulation  is  nearly  identical  to  the 
measured  data.  The  relative  amplitudes  of  the  peaks  in  the  experimental  data 
fall  between  the  results  obtained  with  the  ideal  (flat)  power  spectrum  and 
those  obtaiined  by  weighting  the  simulation  with  the  power  spectrum  of  Figure 
8. 


Figure  30  demonstrates  the  effect  of  temperature  on  the  output  frequency 
of  the  klystron  oscillator.  The  specification  sheet  indicates  that  the  Varian 
VRB-2111A  reflex  klystron  has  a  typical  frequency/temperature  coefficient  of 
-1.6  MHz/°C;  all  of  the  simulation  results  used  this  value.  The  effect  of 
this  coefficient  is  to  introduce  a  constant  phase  shift  into  the  interference 
patterns.  This  mechanism  was  used  to  adjust  the  phase  of  the  simulated  results 
to  agree  with  the  experimental  results,  and  the  temperature  correction  factors 
are  shown  in  the  figures. 

The  phase  of  the  signal  is  really  of  no  importance,  since  the  desired 
Information  is  the  Fourier  transform  of  the  measured  data.  Other  sources  of 
slow  frequency  Instability,  such  as  power  supply  voltages  will  contribute  a 
phase  factor  which  may  be  considered  to  be  constant  for  a  given  measurement. 

The  effect  of  the  nonlinear  frequency  dependence  is  to  cause  a  changing 
periodicity  in  the  interference  patterns.  Figure  13  shows  that  the  slope 
becomes  steeper  on  the  frequency  versus  voltage  curve  at  the  upper  end. 
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INPUT 

f  No.  of  Scatterers 
^  No.  of  Frequency  Points 

a  Estimated  Temperature  Change  (C°) 
"  of  the  Source 


READ-  - 

^  Power  Spectrum  of  the  Source 
(Amplitude  vs  Frequency) 


TNPtTT 

^  RCS  of  each  Scatterer 
$  Range  to  each  Scatterer 


Based  on  the  Power  Soectrum: 

^  The  Correct  Nth  Frequency 

^  The  Nth  Power  Weighting  Coefficient 

A  The  Nth  Response  of  the  Square  Law 
"  Detector  (Sq.  3) 


4  Find  the  Maximum  Square  Law  Response 


f  Normalize  the  N  Responses  by  the 
Maximum 


^  Plot  the  Normalized  Response  versus  Applied 
Voltage  on  the  Source 


Figure  20.  Square-law  detector  response  algorithm 


Figure  21.  Discrete  Fourier  transform  algoritlw. 
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Figure  22  Simulation:  two  reflectors  separated  by  6  meters. 
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Figure  23.  Simulation:  two  reflectors  separated  by  5  meters. 
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Klystron  Reflector  Voltage  Change 
Figure  ?4.  Simulation:  two  reflectors  separated  by  4  meters. 
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Figure  26.  Simulation:  two  reflectors  separated  by  2  meters 
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Figure  28.  Simulation:  three  reflectors  separated  by  4m,  2m, 
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AT  =  11°C 


Klystron  Reflector  Voltage  Change 
Figure  29-  Simulation:  three  reflectors  separated  by  3m,  2m. 
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Klystron  Reflector  Voltage  Change 


Figure  30.  Simulation:  two  reflectors  separated  by  1  meter 
temperature  effects. 
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resulting  in  a  decrease  in  the  period  of  the  oscillation  after  the  reflector 
voltage  change  reaches  60  units*  This  effect  has  been  accurately  modeled  in 
the  simulation  by  plotting  the  results  as  a  function  of  reflector  voltage 
change,  rather  than  as  a  function  of  frequency.  Only  when  the  frequency  curve 
is  perfectly  linear  can  the  experimental  data  be  considered  to  be  an  amplitude 
versus  frequency  plot. 

Figure  31  shows  the  output  of  the  FFT  for  the  data  produced  by  the  simu¬ 
lation  for  a  separation  between  reflectors  of  6  meters.  As  noted  in  the 
figure,  the  result  is  for  a  flat  power  spectrum  and  a  perfectly  linear  sweep 
of  the  source.  Figures  32  and  33  show  the  results  of  the  FFT  of  the  data  pro¬ 
duced  by  the  simulation  for  the  other  cases  of  interest.  All  of  the  plots  are 
normalized  to  the  largest  coefficient  (do  term)  of  the  FFT.  Figures  34 
through  36  display  the  same  cases  for  a  separation  of  4  meters. 

Figures  37  and  38  provide  the  results  of  FFT's  of  the  simulated  data  for 
the  three  reflector  case  without  and  with  degradation.  Figure  38  can  then  be 
compared  to  Figure  39.  which  is  an  FFT  of  the  experimental  data  obtained  by 
manually  digitizing  the  data  and  calculating  the  spectrum.  Figure  40  shows 
the  results  for  experimental  data  for  separations  of  3  meters,  2  meters.  The 
FFT  of  the  actual  data  displays  a  broader  and,  hence,  flatter  spectrum  that 
the  FFT  of  the  simulation  result.  The  combined  effect  of  the  nonlinear  power 
curve  and  nonlinear  frequency  sweep  is  to  spread  the  power  over  the  range 
domain.  These  effects  combine  to  mask  the  expected  range  separation  lines  and 
introduce  others  which  make  the  target  appear  more  complicated. 

VII.  CONCLUSIONS  AND  RECOMMENDATIONS 

For  the  data  collected  with  two  corner  reflectors,  there  is  excellent 
agreement  between  the  measured  range  separations  and  those  calculated  from  the 
radar  results.  There  is,  in  general,  good  qualitative  agreement  between  the 
experimental  interference  patterns  and  those  obtained  from  theoretical  com¬ 
puter  simulation  for  both  two  reflector  and  three  reflector  data. 

Most  of  the  temperature  shifts  required  to  adjust  the  phase  of  the  inter¬ 
ference  pattern  (or,  equivalently,  to  adjust  the  frequency  span  end  points) 
are  not  unreasonable.  The  temperature  in  the  High  Bay  during  the  experiment 
was  low  (~10°C),  while  the  surface  body  temperature  of  the  klystron  during 
normal  operation  is  probably  higher  than  50°C.  In  order  to  conserve  tube 
lifetime,  the  klystron  was  often  turned  off  between  measurements,  occasionally 
for  long  periods  of  time.  Thus,  large  temperature  variations  were  possible. 

In  addition,  some  phase  change  may  be  due  to  regulation  or  inadvertent 
adjustment  of  the  klystron  beam  (resonator)  voltage.  An  increase  in  beam 
voltage  causes  a  decrease  in  frequency;  a  crude  measurement  indicated  a  sen¬ 
sitivity  on  the  order  of  -1  MHz/V. 

Two  factors  contributed  to  the  discrepancies  between  the  theoretical  and 
experimental  interference  patterns  with  respect  to  amplitudes:  (1)  the  inabi¬ 
lity  to  adequately  characterize  the  power  curve  (note:  actually  the  power 
curve  used  should  have  been  the  interference  pattern  obtained  with  only  one 
reflector,  since  receiver  characteristics  are  also  included);  and  (2)  because 
the  antenna  beamwidths  were  so  narrow,  the  unlikelihood  that  the  corner 


reflectors  were  equally  Illuminated  by  the  transmit  antenna  or  equally  sub¬ 
tended  by  the  receive  antenna  (especially  for  the  measurements  with  three 
reflectors).  Discrepancies  with  respect  to  the  positions  of  maxima  and  minima 
along  the  "klystron  reflector  voltage  change"  axis  result  from  the  inability 
to  adequately  obtain  the  klystron  frequency  dependence  on  reflector  voltage; 
spectrum  analysis  is  required. 

The  agreement  between  the  Fourier  transform  experimental  and  theoretical 
data  is  less  apparent,  again,  due  to  the  problems  cited  above.  Such  agreement 
is  critical  if  the  sensor  is  to  be  used  to  investigate  complex  targets.  The 
fine  structure  in  the  experimental  data  show  up  as  spectral  lines  which  indi¬ 
cate  wide  range  separations. 

It  is  evident  that  the  basic  radar  measurement,  i.e.,  the  interference 
pattern,  contains  target  information  greatly  modified  by  the  dynamic  sensor 
characteristics.  In  order  that  the  data  are  more  representative  of  the  target 
than  of  the  radar,  the  following  radar  attributes  are  deemed  desirable: 

1.  The  frequency  of  the  source  should  be  a  linear  function  of  the 
control  voltage,  or  the  control  voltage  be  such  that  the  frequency  sweep  is 
linear,  and  that  a  voltage  proportional  to  frequency  be  provided. 

2.  The  transmitted  and  received  power  curve  should  be  flat  over  the  fre¬ 
quency  diverse  bandwidth,  or  it  should  be  adequately  characterized  and 
allowances  made  for  it  in  signal  processing;  this  requirement  is  increasingly 
more  critical  for  increasingly  more  complex  targets. 

3-  The  frequency  diverse  bandwidth  should  be  as  large  as  possible. 

4.  The  transmitter  and  receiver  should  be  temperature  controlled,  or 
their  temperature  characteristics  should  be  determined  over  the  operating 
ranges  and  compensated  for.  .'ote  that  a  phase  change  in  the  interference  pat¬ 
tern  does  not  affect  the  Fourier  transformed  data,  provided  that  the  bandwidth 
is  sufficiently  large. 

5*  A  pulsed  system  is  preferable  to  a  continuous  wave  system,  in  order 
to  limit  the  effects  of  clutter  (such  as  occurred,  as  a  result  of  the  High  Bay 
door).  Pulse  width  can  be  adjusted  to  maximize  target  response  while  mini¬ 
mizing  clutter  contamination.  Pulse-to-pulse  coherence  is  not  required. 

In  terms  of  bandwidth,  temperature  effects,  power  curve  and  frequency 
nonlinearities,  and  continuous  wave  operation,  the  94  GHz  reflex  klystron  is 
clearly  not  an  adequate  source  for  application  to  complex  targets.  The  radar 
at  140  GHz  had  a  reasonable  bandwidth,  but  also  had  severely  nonlinear  power 
characteristics.  It  is  recommended,  therefore,  that  further  research  be  per¬ 
formed  in  an  effort  to  find  suitable  millimeter  wave  sources  and  to  build  and 
test  experimental  systems  achieving  the  desirable  characteristics  delineated. 

Three  millimeter  wave  (near  3mm)  sources  are  currently  available  which 
have  large  electronically  tuned  frequency  diverse  bandwidths.  The  first,  a  94 
GHz  solid  state  source  employing  a  varactor  controlled  gunn  diode  oscillator, 
has  a  bandwidth  of  600  MHz  with  a  continuous  wave  power  output  of  30mW. 
Although  this  power  level  is  probably  too  low  for  tactical  systems,  the  device 
is  coherent,  and  pulse  compression  may  be  employed.  The  cost  is  about  $16  k. 


The  second  source  is  a  80  GHz  rising  sun  magnetron  tube  having  a  3-4 
kilowatt  peak  power  output  and  a  1  GHz  agile  bandwidth.  Rapid  tuning  is 
achieved  by  movement  of  a  small  conducting  ring  at  one  end  of  a  rising  sun 
resonant  cavity.  The  tuning  ring  is  supported  and  moved  by  a  piezoelectric 
transducer.  The  magnetron  pulse  duration  is  typically  50  ns  with  a  duty  cycle 
of  0.0002.  A  high  voltage  power  supply  (lOOkW)  is  required.  The  cost  of  the 
tube  is  only  about  $20  k. 

The  third  off-the-shelf  source  is  an  0-type  backward-wave  oscillator 
(carcinotron) ,  electronically  tunable  from  88  GHz  to  92  GHz.  Continuous  wave 
power  over  this  bandwidth  is  0*3  to  4  watts,  although  higher  power  can  be 
achieved  with  narrower  bandwidths.  Tube  cost  is  $20-25  k.  Kilovolt  power 
supplies  are  required. 

In  addition,  there  is  in  developement  a  94  GHz  coupled  cavity  traveling 
wave  tube  (TWT)  with  a  2  GHz  bandwidth.  Power  output  is  100W  with  a  0.50 
duty  cycle. 
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2  M  D 1 11  r:  1  i 00  J  j  GC  256  ] »  d.  256  I » NC  256  I 
40  i)TSC  NUMBER  OF  SCATTERS*" 5 
50  INPI.it  N 

60  MSP  "ENTER  #  UF  STEPS"! 

70  INPUT  N 1 

00  DfSp  "ENTER  TEMP  CHANGE < C  3 "  5 

90  INPUT  Ci 

ISO  FOR  1=1  TO  2i 

130  RE+U)  Cl  I J 

140  HE,,,  1 

150  FOR  ;  -■!  TO  21 

160  READ  Pi  i  1 

170  NEXT  1 

ISO  FOP  L-J  TO  N 

190  DlSF  'ErNiER  Res  FUR  EACH  SCAT  TER  ' 5 
200  DI SR  L» 

210  INPUT  SCI.] 

226  Hi  SR  'ENTER  RANGES  Ml  TO  THE  SCHTTERER"! 

230  I: ISP  l 5 
248  INPUT  PEL.] 

256  NEXT  L 

269  FOR  K>  1  TO  HI 

270  2=!'  K~  i  1  •••' 0  H 1  ■■"20 > 

280  Zt~Z-lNT(Z> 

290  IF  21\  IE -06  THEN  330 

300  Gr  K  3=  •"  Ft  I  NT  C  7 3+2  3  -T  [  I  NT  (  23  +  1  J  ?*2 1 +F  [  I  NT  (23  +  13+  <  <  - 1 . 6E+06  )  *C  1 3 
3 1 0  CE  K  3=(  HE  I  NT-.  2  '• +2  3-  AE  I  NT  (2)  +  l  3  >  *2 1 +RE  I  NT  (23  +  13 
320  GOT 0  350 

330  CE  K  3=FE  I  NT  ( 2  >  t  1  3+ <  - 1 . 6E+06+C 1 3 
340  Oi  K  ] --At  TNT <23  +  1  3 


-  V 


3  6  0  FOR  I  - 1  T  0  ('  H  - 1 ) 

378  FOR  J=<I+1>  TO  N 
3S0  R 1 =R[ J  J— R[  I  ] 

390  C=3E+0S 

400  fll=4*PI *GC K ]* ( R 1 /C > 

4 1 0  B- < CC K ] ) * < SC  1 J+St J ]+2*SQR < SC  I ]*SC J ] > *COS <fll> > 

420  V1=V1+B 

430  NEXT  J 

440  NEXT  I 

450  VC K  ]=V 1 

460  IF  K > 1  THEN  490 

470  V3=VC K  ] 

480  GOTO  510 

490  IF  VCKKV3  THEN  510 

5G0  V3=VC K  ] 

510  NEXT  K 
520  SCALE  0  ,  N 1 , 8 , 1 
530  FOR  K= 1  TO  N1 
540  PLOT  K.t  VCK  ]/V3> +0 
550  NEXT  K 
560  S TOP 

570  CAT A  0 1 0  j  O . U 0 . 2 1 j 0 . 375 , 0 . 53 j  0 . 67 , 0 . 785  , 0 . 875  j  0 . 935  > 0 . 98  >  0 . 99  j  8 . 95  , 0 . 88 
580  CAT  A  0 .765*0. 595 >0.4*0. 1 8 > 0 . 8 1 >  0  > 0 

590  CAT A  9 . 3953E+ i 0 > 9 . 39655E+ 1 0 >  9 . 39735E+ 1 0,9. 39825E+ 1 0,9. 399E+ 10,9. 3997E+ 1 0 
60O  CAT A  9 .  4085E+ 10,9.401 3E  + 1 0,9. 402 1 E+ 1 0 >  9. 40285E+ 10,9. 4035E+ 10,9. 40425E  + 1 0 
6 1 0  CRT A  9 . 405E+ 1 0 , 9 . 4859E+1 0,9. 4069E+1 0,9. 4079E+ 10,9. 40905E+ 1 0,9.41 055E+ 1 0 
620  CHTA  9. 412E+10, 9. 41345E+10, 9. 4149E+10 
1500  END 
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ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  r  o  r  o  r  o  ro 


REM  THIS  PROGRAM  WILL  TAKE  A  N  POINT  BFT  OF  AH  INPUT 
REM  THE  DATA  ARRAY  MUST  BE  REAL.  THE  NUMBER  OF  DATA 
REM  OF  2,  NOT  TO  EXCEED  8  <?56  POINTS).  THIS  PROGRAM 
REM  ARRAYS  ARE  DIMENSIONED  IN  THE  CALLING  PROGRAM. 

REM  THIS  IS  AN  IN  PLACE  f FT.  THE  ELEMENTS  OF  DATA  <V< 
REM  BY  THE  FFT  COEFFICIENTS. 

N2=N1.  2 


DATA  ARRAY  V I  > 
POINTS  MUST  A  POWER 
ASSUMES  THAT  ALL 

I>>  ARE  REPLACED 


N3=H 1-1 
J=1 

REM  G( I )=REAL  PART  OF  V(I) 

REM  C< I >  =  IMAGINARY  PART  OF  V<1> 
FOR  1=1  TO  HI 


Cl  I 1=0 
Gt I  1=VC  I  1 
NEXT  I 
GUSUB  2T'Y0 

REM  THE  ABOVE  STATEMENTS  INITIALIZE  YU)  &  2(1  > 
REM  THE  #  OF  POINTS  MUST  SATISFY  N=2fM 
DISP  "INPUT  M "  j 
INPUT  M 
XI  =  0 

FOR  L= 1  TO  M 
L2=  INI  (2t'L) 

L3=  I  NT  ( L2-  2 ) 

U 1  =  1 
U2=0 

U1=C0S(PI/L3> 

W2=SIN-::PI/L3) 

FOR  J=1  TO  L.3 

FOR  I  =  J  TO  HI  STEP  L.2 

1 1  =  INT  < I +L3) 

T 1  -  <GC  II  ]*UlWe[  II  1*U2) 

T2=<G[  II J+U2 )  +  ( CC II  ]*U1> 

GI  Ii ]=GC I  1-T 1 
CC II ]*CC I  1-T2 
GC  I  1=GI  I  I+T 1 
CL  I 1=CC I  1+T2 
NEXT  i 
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2371  H=U1 

2372  31=U 2 

2  3 8 0  U 1  ( n  *  W 1  >  -  ( Ii  1*W2  > 

2  3 9 3  IJ 2  =  <  D * W 2 )  +  ■:  D  1*W1  > 

2489  NEXT  J 
2410  NEXT  L 

2420  FOR  1=1  TO  N1 

2430  XC I 3=SGR ( (GC I 3*2>+<C[ I  ]t 

2440  IF  XCIKX1  THEN  2460 

2450  Xl=XtII 

2460  NEXT  I 

2480  SCALE  0>N1>0j 1 

2490  FOR  1=1  TO  Hl+1 

2500  PLOT  1-1 j  Xt I  3  "XI j  +2 

2501  I  PLOT  1 j -XC I  ]/X 1 j  + 1 
2510  NEXT  I 

2520  STOP 

2570  FOR  1=1  TO  N3 
2580  IF  I  >=  J  THEN  2620 
2590  T 1 =G[ J  ] 

2600  GC J  ]=G[  I  ] 

2610  Gt I  ]=T 1 
2620  K=N2 

2630  IF  K  >=  J  THEN  2678 

2640  J-J-K 

2650  K=K  '2 

2660  GOTO  2630 

2670  J= J+K 

2680  NEXT  I 

2690  RETURN 

2700  END 


10  H I M  VI  256  3*  SC  25  3 ,  RC  25  3 » fl[  100] 

20  D I M  FT  100  Is  GC  256  3 »  CC  256  ] ,  XC  256  3 
40  LISP  "HUMBER  OF  SCATTERS*" 5 
50  INPUT  N 

80  DISP  "ENTER  TEMP  CHANGE ( C > " 5 

90  INPUT  Cl 

120  FOR  1=1  TO  21 

130  READ  At.' I  3 

140  NEXT  I 

150  FOR  1=1  TO  21 

160  READ  FC  13 

170  NEXT  I 

171  DISP  "ENTER  #  OF  STEF'S"  5 

172  INPUT  HI 

180  FOR  L= 1  TO  N 

196  DISP  "ENTER  RCS  FOR  EACH  SCATTER  '  5 
200  DISP  Li 
210  INPUT  SC  L  3 

220  DISP  "ENTER  RANGE CM)  TO  THE  SCATTERER" i 

230  DISP  Li 

240  INPUT  RC  L  3 

250  NEXT  L 

260  FOR  K=l  TO  HI 

270  Z=CK- 1 ) / <. N 1  /20 > 

280  21-Z-INT C2) 

290  IF  Z1C1E-06  THEN  330 

300  GC  K 3=  C  FC I  NT  < Z  >  +2 3-FC I NT  C Z  >  + 1 3 ) *Z 1 +FC I NT C  Z  >  + 1 3+  C  < - 1 . 6E+Q6 ) *C 1 
310  CC  K 3= (AC  I  NT  CZ)+2  3- AC INTCZ)  +  1 ]>*Z1+AC I NT  <  2  >  + 1 3 
320  GOTO  350 

330  UC  k  3=FC  I  NT  (.  Z )  + 1  3+ C  ~  1 . 6E+06+C 1 ) 

340  CC  K.  3= AC  I  NT  CZ)  +  1  3 
350  V1=0 
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FOR  1  =  1  TO  < N  -  1 ) 

FOR  J=(I+1>  TO  N 
R 1  =R[  J ]-R[  I  I 
O3E+08 

R 1 — 4  *  P I  G  C  K  ]*(.  R1."C> 

o-a;c k  i>*< sc  i  j+sr  j ]+2*sqr<'SC  i  ]*sc  j  :i>*cos<ai  > > 

V 1 =V 1 +B 
NEXT  J 
NEXT  I 
VLR  ]=V1 

IF  K > 1  THEN  490 
V3=VC K  ] 

GOTO  510 

IF  VCKKV3  THEN  51 0 
V3-VC K  ] 

NEXT  K 

GOTO  560 

SCALE  0  ?  N 1  ?  0 j 1 

FOR  K=1  TO  Ml 

PLOT  K?VCK]/V3>+0 

NEXT  K 

STOP 

DAT A  8 ?  0  j  0 , 1 ? 0 . 2 1 j  0 . 075  * G . 53  ? 0 . 67? 0 . 785  ?  0 . 875 j  0 . 935 ? 0 . 98  ? 0 . 99  ? 0 . 95 ? 0 . 88 
DATA  0 . 765 ? 0 . 595  >0.4?0. 1 8  ?  0 . 0 1?0?0 

DATA  9. 3958E+10? 9. 39655E+10? 9. 39735E+10? 9. 39825E+10?  9. 399E+10?  9. 3997E+10 

DAT A  9. 4005E+1 0? 9. 4013E+10? 9. 4021E+10? 9. 40285E+10? 9. 4035E+10? 9. 40425E+1O 

DATA  9 , 405E+10? 9. 4059E+10? 9. 4069E+10? 9. 4079E+10? 9. 40905E+10? 9.41055E+10 

DA T A  9.41 2E+ 10?  9. 4 1 345E+ 10?  9 . 41 49E+ 1 0 

LINK  2? 640 

N2=N  1  •'  2 

N3=N 1  - 1 

J=1 

REM  Gil) "REAL  PART  OF  VCD 

REM  C(I)  =  IMAGINARY  PART  OF  VCD 

FOR  1=1  TO  HI 

C[ I  ]=0 

bC I  3=VE  I  I 

NEXT  I 

G OS UB  1160 

REM  THE  ABOVE  STATEMENTS  INITIALIZE  YCI>  &  Z(I) 

REM  THE  #  OF  POINTS  MUST  SATISFY  N=2TM 
DISP  "INPUT  M"? 

INPUT  M 
X1=0 

FOR  L=1  TO  M 
L2=INT  <2t'L> 
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